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We study the breakup of one-neutron halo nuclei in the Coulomb field of a 
target nucleus. In the post-form distorted wave Born approximation theory of 
this reaction, with only Coulomb distortions in the entrance and the outgoing 
channels, an analytic solution for the breakup T-matrix is known. We study 
this T-matrix and the corresponding cross-sections numerically. This formula 
can be related to the first order semiclassical treatment of the electromagnetic 
dissociation. This theory contains the electromagnetic interaction between the 
core and the target nucleus to all orders. We show that higher order effects 
(including postacceleration) are small in the case of higher beam energies 
and forward scattering. We investigate the beam energy dependence of the 
postacceleration effects. They are found to be quite important for smaller 
beam energies (slow collisions), but almost negligible at larger ones. 
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I. INTRODUCTION 



Breakup processes, in nucleus-nucleus collisions, are complicated, in whatever way they 
are studied. Coulomb breakup (CB) is a significant reaction channel in the scattering of 
halo nuclei from a heavy target nucleus (see, e.g., With the operation of exotic beam 

facilities all over the world, these reactions (previously restricted essentially to deuteron 
induced reactions) have come into sharp focus again. CB provides a convenient way to 
put constraints on the structure of these nuclei This is of interest also for nuclear 

astrophysics, since the breakup cross section can be related to the photo-dissociation cross 
section and to radiative capture reactions relevant for nuclear astrophysics [0]. 

The breakup reactions of the halo nuclei have been investigated theoretically by several 
authors, using a number of different approaches (see, e.g., for an extensive list of 

references). One often used method has been to treat the breakup reaction as the inelastic 
excitation of the projectile from its ground state to the continuum [^. The corresponding 
T-matrix is written in terms of the prior form distorted wave Born approximation (DWBA) 
||1CI||. For pure Coulomb breakup, the semiclassical approximation of this theory is the first 



order perturbative Alder- Winther theory of Coulomb excitation W^. It has recently been 



used to analyze the data on the breakup reactions induced by the one-neutron halo nuclei 
^^Be [O and ^^C [O]. The question of the higher order electromagnetic effects fT^-ITH 



was 



studied recently within this framework. They were found to be small, for both zero 
range as well as finite range wave functions of the relative motion of the fragments within 
the projectile. Also, in a dynamical description of the breakup of these nuclei where the time 
evolution of the projectile in coordinate space is described by solving the three-dimensional 
time dependent Schrodinger equation treating the projectile-target interaction as a time 
dependent external perturbation, the higher order effects turned out ||20| to be only of the 
order of 10% for beam energies in the range of 60 - 80 MeV/nucleon. 

A direct breakup model (DBM) (which reduces to the Serber model in a particular 
limit 0) has been formulated within the framework of the post form DWBA An 
important advantage of this model is that it can be solved analytically for the case of 
the breakup of the neutron halo nuclei with the entrance and outgoing channels involving 
only the Coulomb distortions It constitutes an ideal "theoretical laboratory" to 

investigate the physics of the breakup reactions, its certain limiting cases, and its relation to 
other models like the semiclassical approximation. Particularly, the effect of postacceleration 
(to be explained in more detail below) can be studied in a unique way within this approach. 

The aim of this paper is to investigate the role of the higher order effects (which includes 
postacceleration) in the Coulomb breakup of the one-neutron halo nuclei, ^^Be and ^^C, 
within the post form DWBA theory of the breakup reactions. We calculate the triple and 
the double differential cross sections of the fragments, observed in the breakup of these nuclei 
on a 208p]^ target, within the exact theory as well as within its first order approximation. We 
also calculate the relative energy distributions of the outgoing fragments emitted in these 
reactions, within the two theories. Calculations have been performed for a range of beam 
energies in order to investigate the beam energy dependence of the higher order effects. 
In Sec. II, we present our theoretical model. The results of our calculations and their 
discussions are given in Sec. III. Summary, conclusions, and the outlook of our work in 
given in Sec. IV. 
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II. DESCRIPTION OF THE THEORETICAL MODEL 



A. Coulomb wave Born approximation 

We consider the reaction a-\-t — > c + n + t, where the projectile a breaks up into the 
charged fragment c and the neutron n in the Coulomb field of a target t. 

The starting point of the Coulomb wave Born approximation (CWBA) is the post form 
T matrix of the breakup reaction which is given by 

The functions xi~^*(kc, r), and x^r^*(k„, r„) are the distorted waves for the relative motions 
of c and n with respect to t and the center of mass (cm.) of the c-\-t system, respectively. The 
arguments of these functions contain the corresponding Jacobi momenta and coordinates. 
$'s are the internal state wave functions of the concerned particles which depend on the 
internal coordinates ^. The function ri, r^) is the exact three-body scattering wave 

function of the projectile with a wave vector k^, which satisfies outgoing wave boundary 
conditions. The vectors kc and k„ are the Jacobi wave vectors of c and n, respectively, in 
the final channel of the reaction. The function V^n(ri) represents the interaction between c 
and n. For the pure Coulomb breakup case, the function Xc ''(kc, r) is taken as the Coulomb 
distorted wave (for a point Coulomb interaction between the charged core c and the target) 
satisfying incoming wave boundary conditions, and the function Xn ''(knji'n) is just a plane 
wave as there is no Coulomb interaction between the target and the neutron. The position 
vectors satisfy the following relations 

r = Tj - ari, a = — , (2) 
rria 

r„ = 7ri + (^ri, 5= — ^ — , 7 = (1 - q;5) . (3) 

In the distorted wave Born approximation, we write 

^(+)(ea,ri,r,) ^ $,(ea,ri)xi+nka,r,). (4) 

In Eq. (4), the dependence of on ri describes the relative motion of the fragments c and 
n in the ground state of the projectile. The function xi^^(ka,rj) is the Coulomb distorted 
scattering wave describing the relative motion of the cm. of the projectile with respect to 
the target, satisfying outgoing wave boundary conditions. 

The integration over the internal coordinates ^, in Eq. (1), gives 



/ 



d^K{QK{^n)^a{U ri) = {imjnf^n\jfJ'){jcf^cjf^\jaf^a)i^^a{ri), (5) 



with 

$a(ri) = ue{n)Ye^{h). (6) 
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In Eq. (6), £ (the orbital angular momentum for the relative motion between fragments c 
and n) is coupled to the spin of n and the resultant channel spin j is coupled to spin of 
the core c to yield the spin of a {ja)- The T-matrix can now be written as 

^ ^ {^mJr,^ln\J^^){Jc^^c]^^\Ja^^a)^'lP?J^^^{K, K\ kj, (7) 

where 

^>£™'^(kc, k„; k„) = j drldr,x^)*(ke, v)e-^'^-^-VUn)u,{n)Y,^{h)x^+\K. r.), (8) 
with i = VWTT. 

Eq. (8) involves a six dimensional integral which makes its computation quite compli- 
cated. The problem gets further acute because the integrand has a product of three scat- 
tering waves that exhibit an oscillatory behavior asymptotically. Therefore, approximate 
methods have been used, such as the zero range approximation (ZRA) in which we write, 
Vcn(ri)^'a(ri) = -Do^(ri) with Dq being the zero range normalization constant, (see, e.g.. 



^), or the approximation used in [^, where the projectile cm. coordinate is replaced 
by that of the core-target system (i.e. «i r). Both these methods lead to a factorization of 
the amplitude [Eq. (8)] into two independent parts, which reduces the computational com- 
plexity to a great extent. However, the application of both these methods to the reactions 
of halo nuclei is questionable (see, e.g., 0, for a detailed discussion). 

In the finite range CWBA theory , the Coulomb distorted wave of particle c, in the 
final channel, is written as 

xi-\Kr) = e-''^'^-^^xi-\K,ri). (9) 

Eq. (9) represents an exact Taylor series expansion about r j if K = —iVr, is treated exactly. 
However, instead of doing this we employ a local momentum approximation p6| , p7| , where 
the magnitude of momentum K is taken to be 



K{R) = ^'^{E-V{R)). (10) 

Here m is the reduced mass of the c — t system, E is the energy of particle c relative to 
the target in the cm. system and V{R) is the Coulomb potential between c and the target 
t separated by the distance R. Thus, the magnitude of the momentum K is evaluated at 
some separation R which is held fixed for all the values of r. For further details and the 
discussion on the validity of this approximation, we refer to |8|j2^. 

On substituting Eq. (9) into Eq. (8), we obtain the following factorized form of the 

amplitude (3^^^ 

ippr^^^K K; K) = Zem J dvair>iK r,)e-''''-'^^xi^\K, r,). (11) 



where 



-'£m 



J drie-''''-'^VUri)ue{ri)YeUri), (12) 
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where ki = 7k„ — aK. 

This amphtude differs from that of the ZRA studied earher |]lOl as it allows the use of 
the full wave function for the relative motion of the fragments (corresponding to any value 
of i) in the ground state of the projectile. However, it should be stressed that as for as the 
postaccetelation effects are concerned, both the amplitudes would lead to identical results. 

The triple differential cross section of the reaction is given by 



27T 



dEcdVtcdVtn 



,CWBA|2 



(13) 



where p{Ec,, fie, ^n) is the appropriate three-body phase space factor [§], and Va the velocity 
of particle a. 

On substituting the Coulomb distorted waves. 



X 



(-)*(k„ r,) = e--^^I^T{l + zryJe-^'^-'-»iFi(-z77„ 1, i{Kr, + k,.r,)) 



(14) 



(k„ r,) = e--'T^'''V{l + ir^ay''-'\F^{-ir,a, 1, i^Kn - k^.r,)) , (15) 
into Eqs. (11) and (13), one gets for the triple differential cross section 



laic I ^ 1 2 I 7' 1 2 

dEjnjQn hva (e2-'?= - l)(e2-'?« - 1) ' ' ^ ' ^' ■ 



-p{Ec, fic, fin) 



(16) 



In Eqs. (14-16), rj^s are the Coulomb parameters for the respective particles. In Eq. (16), / 

can be evaluated in the closed foi 

VaVc)2Fiil - iva,l - tVc;^; D{0)) 



is the Bremsstrahlung integral |29] which can be evaluated in the closed form: 

/ = _45(o)f^^ 



\dAJ A=o 
\dAJ A=o 



where 



5(A) 



{k^ - 2k.ka - 2AA;„)^''°(P - 2k.ke - 2Ak, 



DiA) 



2P{kak, + kg.k,) - Ajk.kg + Afc„)(k.k, + Afcj 

(F - 2k.k„ - 2Aka){k^ - 2k.kc - 2AA;e) 



with 



k — kfl kc Skji- 

The factor Z'l contains the projectile structure information and is given by 

Z'e = / drirlje{kiri)Vcn{ri)ui{ri). 



(17) 



(18) 



(19) 



(20) 



(21) 
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It may be noted that the triple differential cross sections with respect to relative and 
cm. coordinates of the fragments are related to those given by Eq. (13) as 

^± = _L21^ (22) 

where the form of Jacobian J is the same as that given in PDj. In Eq. (14), t — {c + n) 
corresponds to the coordinates of the relative motion of the cm. of the fragments c and n 
with respect to the target, while c — n corresponds to that of the relative motion between 
them. 

The CWBA description [Eqs. (11-12)] includes the effects of postacceleration, which 
refers to the situation where the core c has a larger final state energy than what one gets 
from sharing the kinetic energy among the fragments according to their mass ratio. This 
effect arises in a purely classical picture |^ of the breakup process. The nucleus a = (c + n) 
moves up the Coulomb potential, loosing the appropriate amount of kinetic energy. At an 
assumed "breakup point" , this kinetic energy (minus the binding energy) is supposed to be 
shared among the fragments according to their mass ratio (assuming that the velocities of c 
and n are equal). Running down the Coulomb barrier, the charged particle c alone (and not 
the neutron) gains back the Coulomb energy, resulting in its postacceleration. Of course this 
picture is based on the purely classical interpretation of this process, and will be modified 
in a quantal treatment, where such a "breakup point" does not exist. Postacceleration is 
clearly observed in the low energy deuteron breakup, both in the theoretical calculations as 
well as in the corresponding experiments (see, e.g., [0,^). However, in the description of 
the Coulomb dissociation of halo nuclei at high beam energies within this theory [p|, p2| , |33| , 
the postacceleration effects become negligibly small. We shall investigate this point further 
for the ^^Be and ^^C Coulomb dissociation experiments |T^|T^. 

On the other hand, in the first-order semiclassical Coulomb excitation theory which was 
widely applied in the recent years to the Coulomb dissociation of the neutron halo nuclei 
(see, e.g., [0), these effects were found to be small, for both zero range as well as finite 
range wave functions of the c + n system. 

One can establish a relation between the apparently very different CWBA and the semi- 
classical theory. It was recently noticed that in the limit of Coulomb parameter 77^ <C 1 
(i.e. in the Born approximation), both theories give the same result. It was further found 
that this agreement is also valid for arbitrary values of rja and rjc, provided the beam energies 
are high as compared to the relative energy [Ecn) of fragments c and n in the ground state 
of the projectile. The first order approximation to the amplitude given by Eq. (11), can be 
written as 



ofirst order 



4^^^/couie 



X 



rrir 



[kl - {6kn - kaf] 



(23) 



where the relative phase (j) = cr('7c) — o'iVa) — — ^ / [2 log |D(0)|], with a{ri) being the usual 
Coulomb phase shifts, and ^ = rjc — rja. In Eq. (23), we have defined f^Qy^i = ^rjaka/k"^. This 
term is very similar to the Born approximation (BA) result given in in the limit ^ ^ 
it actually coincides with the BA expression. This equation can be used to investigate the 
role of higher order effects. 
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B. Scaling Properties 



In many experimental situations, the momentum transfer k [Eq. (20)] is small. One can 
expand Eq. (23) (with = ^ = 0) for small values of k to obtain 

«^first order _ . 2Z^ 2q ■ k 

^^"^ ~ •'coul ^2 1^ , _ ^3^,^2 I ^2^2' ^ ' 



where the relative momentum between c and n is sriven by o = mnkc ^ jg j-gi^^g^;! 

the c—n separation energy in the ground state of the projectile E^^^'^(= ^^fi;^/2/i, with /i being 
the reduced mass of the c-n system) . This result is in remarkable agreement with the usual 
first order treatment of the electromagnetic excitation in the semiclassical approximation. 

In the semiclassical approach, the scattering amplitude is given by the elastic scattering 
(Rutherford) amplitude times an excitation amplitude a(&), where the impact parameter h 
is related to the transverse momentum transfer, g^, and rjahy h = 2r]ah/q±. The absolute 
square of a{b) gives the breakup probability P{b), which, in the lowest order (LO), is given 

by Hjn 



X 



4 



dq 37rK(l + x2)4^ 



(25) 



where variable x is related to the relative momentum between n and c by x = and y is a 
strength parameter given by 

2Z+Zcm„e^ 

y = . / 1 u • 26 

This formula shows very interesting scaling properties: very many experiments, for neutron 
halo nuclei with different binding energy, beam energy, and scattering angles (or and 
all lie on the same universal curve\ The deviations from this simple scaling behavior, e.g., 
postacceleration effects, will lead to violations of such scahng. 



III. NUMERICAL RESULTS 

We now investigate the breakup of the one-neutron halo nuclei ^^Be and ^^C. We take a 
heavy target of atomic number Z=82. In this paper, all the higher order results correspond 
to calculations performed within the finite range CWBA model, while the first order results 
have been obtained by using Eq. (23). The structure term Z'^ [Eq. (21)] was calculated by 
adopting a single particle potential model to obtain the ground state wave function of the 
projectile. The ground state of ^^Be was assumed to have a 2si/2 valence neutron coupled 
to the ^°Be(0"'") core with a binding energy of 504 keV. The corresponding single particle 
wave function was constructed by assuming the neutron - ^°Be interaction of the Woods 
- Saxon type. For a set of values of the radius and diffuseness parameters (same as those 
given in |^), the depth of this potential was searched so as to reproduce the ground state 
binding energy. This 2si/2 wave function has an additional node as compared to a simpler 
zero-range wave function. For the ^^C, the ground state was assumed to have a [^^C(0"*") 
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®2si/2^] configuration with a separation energy of 0.530 MeV. Tlie radius and diffuseness 
parameters, used in tlie well-deptli searcli, were tlie same as tliose given in 

In Fig. 1, we present calculations for the triple differential cross sections for the breakup 
reaction ^^Be + Pb ^ n + ^"Be + Pb, as a function of the energy of the ^"Be core (E^), 
for four beam energies lying in the range of 5 MeV/nucleon - 72 MeV/nucleon. To see the 
postacceleration in a clear way, it is very useful to study the cross-section as a function of 
the core energy. The results obtained within the higher order and the first order theories 
are shown by solid and dotted lines, respectively. 

It can be seen from this figure that while for lower beam energies, the higher order and 
first-order results differ considerably from each other, they are almost the same for the beam 
energy of 72 MeV/nucleon. In each case, the first order cross sections peak at the energy of 
the core fragment which corresponds to the beam velocity (this value of the core fragment 
energy will be referred to as E;,^ in the following). In contrast to this, the peaks of the higher 
order cross sections are shifted to energies > E^^ for the three lower energies. Only for the 
72 MeV/nucleon beam energy, does the higher order result peak at Ef,^. This shows very 
clearly that the finite range CWBA model exhibits postacceleration for beam energies < 
30 MeV/nucleon, while this effect is not present at 72 MeV/nucleon. Therefore, the higher 
order effects are minimal for the Coulomb breakup of ^^Be at the beam energies > 70 MeV. 



This result is in agreement with those obtained in 1 19,20 



> 

CD 



E 
o 
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FIG. 1. Triple differential cross section as a function of the energy of Be core for the reaction 
"Be + Pb ^ n + ^^Be + Pb at the beam energies of 72 MeV/nucleon, 30 MeV/nucleon, 10 
MeV/nucleon and 5 MeV/nucleon. The results of the finite range CWBA and first-order theory 
are shown by solid and dotted lines respectively. 
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In Fig. 2, we compare the first-order (dotted lines) and the higher order (sohd hnes) 
resuhs for the double differential cross section for the same reaction, and for the same beam 
energies as in Fig. 1. These results have been obtained by integrating the triple differential 
cross sections over the unobserved neutron sohd angles. The cross sections are shown as a 
function of Ep. It is clear that for beam energies < 30 MeV/nucleon, the first order results 
peak at Ef,^, hut the higher order cross sections have their maxima at energies larger than 
Eb^. In contrast to this, both the higher order and the first order cross sections peak at 
the same value of Ec (= E^^), for the 72 MeV/nucleon case. Therefore, the postacceleration 
effects are pronounced for the smaller beam energies, whereas they become quite small for 
the higher energies. The near equality of the first-order and the finite range CWBA cross 
sections, at the beam energy of 72 MeV/nucleon, suggest that for this reaction, the higher 
order effects, in general, are quite irrelevant at beam energies > 70 MeV/nucleon. 
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FIG. 2. Double differential cross section as a function of the energy of ^°Be core for the 
reaction ^^Bc + Pb ^ n + ^^Bc + Pb at the beam energies of 72 McV/nucleon, 30 MeV/nucleon, 
10 MeV/nucleon and 5 MeV/nucleon. The results of the finite range CWBA and first-order theory 
are shown by solid and dotted lines respectively. 

In Fig. 3, we compare the results of the first-order and the finite range CWBA calculations 
for the relative energy spectrum of the fragments emitted in the breakup reaction of ^^Be 
on a ^°^Pb target for the same four beam energies as shown in Fig. 2. These cross sections 
have been obtained by integrating over all the allowed values of the angles ^c-n- In both the 
models, the integrations over have been carried out between 1° to grazing angle, 

in the upper two figures, and between 5° to grazing angles, in the lower two figures. The 
integrations over <^t-{c+n) angles have been done over all of its kinematically allowed values. 
The dotted and solid lines represent the results of the first-order and the higher order 
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FIG. 3. The differential cross section as a function of the relative energy of the fragments 
(neutron and ^'^Be) emitted in the ^^Be induced breakup reaction on a ^'^^Pb target at the beam 
energies of 72 MeV/nucleon, 30 MeV/nucleon, 10 MeV/nucleon, and 5 MeV/nucleon. The dotted 
and full lines represent the first-order and the finite range CWBA results, respectively. 

calculations, respectively. We notice that while for the beam energy of 72 MeV/nucleon, the 
higher order effects are minimal, they are quite strong for the lower beam energies, being 
largest at the beam energy of 5 MeV/nucleon. This reinforces the point, already made in 
19| , |20[| , that at the beam energy of 72 MeV/nucleon, the higher order effects are quite small 



if both the first order and the higher order terms are calculated within the same theory. 

In Fig. 4, we show the same results as in Fig. 3, but for the ^^C induced reaction 
on the 208p]^ target for the beam energies of 67 MeV/nucleon, 30 MeV/nucleon and 10 
MeV/nucleon. Like Fig. 4, the integrations over the fragments center of mass angles is done 
in the range of 1° to grazing angle, for first two beam energies and between 5° to grazing 
angle, at the lowest beam energy. We see that in this case too the higher order effects 
are quite weak for the beam energy 67 MeV/nucleon, but appreciable for the lower beam 
energies. 

It may be noted that by comparing the result of a conceptually different model of the 
Coulomb breakup reactions [BH] than ours, with that of the first order semiclassical pertur- 



bation theory of the Coulomb excitation, it has been concluded in |]T8[ that the higher order 
effects are substantial for these reactions even at the beam energies ~ 70 MeV/nucleon. 
However, one should be careful in drawing definite conclusions about the role of the higher 
order effects from such an approach. For a reliable assessment of the contributions of the 
higher order effects, it is essential that both the first order and the higher order terms should 
be calculated within the same model, as has been done in this work (and also in [[191 , ^ ). 
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FIG. 4. The differential cross section as a function of the relative energy of the fragments 
(neutron and ^^C) emitted in the ^^C induced breakup reaction on a ^'^^Pb target at the beam 
energies of 67 MeV/nucleon, 30 MeV/nucleon, and 10 MeV/nucleon. The dotted and full lines 
represent the first-order and the finite range CWBA results, respectively. 



IV. CONCLUSION AND OUTLOOK 

In this paper, we investigated the breakup of the one-neutron halo nuclei, ^^Be and ^^C, 
in the Coulomb field of a heavy target nucleus, within a theory which is formulated in the 
framework of the post form distorted wave Born approximation. This theory contains the 
electromagnetic interaction between the core and the target nucleus to all orders. An attrac- 
tive feature of this formulation is that the corresponding pure Coulomb breakup amplitude 
can be expressed in an analytic form. We have also extracted the first order approximation 
of the full pure Coulomb breakup T-matrix. This enables us to investigate the respective 
roles of the first order and the higher order effects within the same theoretical model. We 
studied the beam energy dependence of the first order and higher order triple and double 
differential cross sections. We also investigated the contributions of the higher order effects 
to the relative energy distribution of the fragments. 

In the higher order model, the peaks in the triple and double differential cross sections 
vs. core energy spectra, are shifted to energies larger than those corresponding to the 
beam velocity, at the incident energies < 30 MeV/nucleon. Therefore, postacceleration 
effects are important at these beam energies. On the other hand, at the beam energy ~ 
70 MeV/nucleon, the corresponding spectra peak at the beam velocity energies, which is 
consistent with no postacceleration. In contrast to this, the first-order cross sections always 
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peak at the beam velocity energy, which is expected as the postacceleration is a higher order 
effect. 

The higher order effects are also found to be quite important in the relative energy 
spectrum of the fragments at beam energies < 30 MeV/nucleon, while they are insignificant 
at the beam energies ~ 70 MeV/nucleon. This suggests that the conclusions arrived at in 
Refs. ||T2| , p^ , where the data on the relative energy spectra of the fragments taken in the 
breakup of ^^Be and ^^C at beam energies ~ 70 MeV, have been analyzed within the first 
order theory of the Coulomb excitation, may not be affected by the higher order effects. 

The present model can be seen as a "theoretical laboratory", which allows to study 
numerically the relation between quantal and semiclassical theories, and the importance of 
postacceleration effects. It should be noted that from an experimental point of view, the 
postacceleration effects are not fully clarified (see, e.g., [p!2| , p7| , |38| ) . Finally, let us mention the 



recent work on the electromagnetic dissociation of unstable neutron-rich oxygen isotopes [|39 



These authors deduce photoneutron cross-sections from their dissociation measurements. If 
the neutrons are emitted in a slow evaporation process in a later stage of the reaction, the 
question of postacceleration is not there. On the other hand, for the light nuclei there is 
some direct neutron emission component and the present kind of theoretical analysis further 
proves the validity of the semiclassical approach used in [P^j. 

Postacceleration effects are also of importance for the use of Coulomb dissociation for the 
study of radiative capture reactions of astrophysical interest. We expect that our present 
investigations will shed light on questions of postacceleration and higher order effects in 
these cases also. 
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